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We have analysed the axisymmetric and non-axisymmetric modes of a con- 
tinuum of vortices in a rotating superfluid. We have investigated how chang- 
ing the temperature affects the growth rate of the disturbances. We find 
that, in the long axial wavelength limit the condition q = a/(l — a') = 1, 
Q\ • where a and a' are temperature- dependent mutual friction parameters, is the 

crossover between damped and propagating Kelvin waves. Thus at temper- 
s' ', atures for which q > 1, perturbations on the vortices are unlikely to cause 
<^> ' vortex reconnections and turbulence. These results are in agreement with 

the recent discovery of Finne et afi* of an intrinsic condition for the onset 
of quantum turbulence in 3 He-B. 



PACS numbers: 67.40. Vs, 67.57.-z 



T3. 

1. MODEL 

O 

When superfluid helium is rotated, an array of quantised vortex lines 

appear which are aligned parallel to the axis of rotation. This array exhibits 
/\ ' oscillation modes which were first predicted and observed by HalP. The 

j^ ■ modes are similar to classical inertial waves in the long wavelength limits and 

are related to the waves of isolated vortex lines in the short wavelength limit. 
These vortex waves are a very important phenomena in the understanding 
of quantised vortex lines and for a review of work in this area see Donnelly 3 . 
The aim of this paper is to extend previous worlP' to consider the tor- 
sional oscillations of a rotating column of superfluid at a range of temper- 
atures. We consider 3 He-B rotating with angular velocity Q = £le z in a 
cylinder of radius r = a. The basic state will consist of a uniform vortex 
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lattice with a large density of vortices aligned along the direction of rotation. 
In this context the behaviour of the superffuid is described by the Hall-Vinen 
equations^. Thus, rather than individual vortex lines, we consider a contin- 
uum of vortex lines. The advantage of this model is that it allows us to 
explore effects which the theory of a single vortex filament 6 cannot describe, 
notably the presence of boundaries and the degrees of freedom represented 
by the coherent oscillatory motion of many vortices. 

The equations of motion of the superfluid in a coordinate system rotat- 
ing with angular velocity fi = $le z may be written as 

£) S 

-J- + (V • V)v s =V$+2u s xn + aAx[Ax (v s - v n )] 

dt „ (1) 

+a'X x (v s - v n ) - av s X x (A-V)A + u s (l - c/)(A-V)A 

where v s and v n are the superfluid and normal fluid velocities in the rotating 
frame, A = V x v s + 2$7, A = A/|A| is the unit vector in the direction of 
A and ^ is a collection of scalar terms. Given the high viscosity of 3 He-B, 
we assume that the normal fluid is in solid body rotation around the z-axis, 
thus in the rotating frame v n = 0. Equation (J2) must be solved under 
the condition that V-u s = 0. The quantity v s = (T/47I") log(6 /a ) is the 
vortex tension parameter, V is the quantum of circulation, a Q is the vortex 
core radius and b = (|A|/r) -1 ' 2 is the average distance between vortices. 
The unperturbed vortex lattice corresponds to the basic state Vq = Vq = 0, 
V^o = for which A = 2£le z , working in cylindrical coordinates (r, (j), z). 

We perturb the basic state by letting v s = u = (u r , u<j,, u z ), 9 = ^o + 
ip, where \u\ <^. 1 and \ip\ <C 1 and linearise the resulting equations of 
motion. We assume normal modes for the perturbations of the form (u, ip) = 
(u(r), ip(r)) exp(iat+im4>+ikz) where m and k are respectively the azimuthal 
and axial wavenumbers and a is the growth rate. The aim of our calculation 
is to determine the real and imaginary parts of a, namely TZe(a) and Tm(o). 

The solution for u z which is regular as r — ► is the Bessel function of 
the first kind of order m, J m ((3r), where 

2 = -/c 2 [(i ( T + ar ? ) 2 + (l-aQV] 

P ~ [(m + ar])(ia + k 2 u s a) + (I - a') 2 u s k 2 ri} {) 

and r\ = 20 + v s k 2 . To determine (3 we enforce the boundary condition 
u r = on the wall of the container r = a, which yields the secular equation 



(ka) 2 J^((3a) 
((3a) J m (f3a) 



(ict + or?) 2 + (1 - a'Yv] + 2mQa{l - a') = 0, (3) 



where the prime denotes the derivative of the Bessel function with respect 
to its argument. 
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Fig. 1. Plots of (a) the growth rate, —Zm{a) (b) \Im(a)/TZe(a)\ of the 
axisymmetric mode (m = 0) against ak for various values of & at T/T c = 0.5. 
The arrows show the direction of increasing £,-, where j = 0, 1, 2, 5, 10, 20, oo. 



2. RESULTS 

Equations (|2|3[) simplify in the axisymmetric case, m = 0. From © we 
find that (5 must satisfy Jq(Po) = Ji(P a ) = 0- Taking £j to be the j th zero of 
J^) = J*(£) (that is Co = 0, Ci = 3.83171, 6 = 7.01559 etc), Equation © 
becomes 



£ 2 



[(ia + ar?) 2 + (1 -a') 2 rj 



2^21 



(fca) 



[(ic + 077) (i<7 + av s k 2 ) 



1 



a 



A2 



ZV^A; 2 ] 



(4) 



This quadratic equation may be solved analytically^! and has an infinite 
number of solutions, according to the value of £j considered. All the solu- 
tions are such that Im(a) is non-negative, so the system is always stable to 
infinitesimal disturbances. The least stable mode is the one for which Im{a) 
is minimum and is given by 



a = ia(fi + v s k 2 ) ± y/{(l - a f ) 2 u s k 2 (2ft + u s k 2 ) - a 2 n 2 } (5) 

which corresponds to £j — > oo. 

Plots of the growth rate, —Im(a) and the ratio of the decay rate of 
the wave to its angular frequency, \Xm(a) /TZe(a)\ can be found in Fig. ^ 
where the mutual friction parameters correspond to T = 0.5T C for which 
a = 0.38 and a' = 0.28, where T c is the critical temperature and we have 
taken v s /a 2 = 0.002 and = 1. The kinks in the curves of Fig. ^a) for 
j > correspond to a switching from being purely imaginary (for which the 
decay of the mode is monotonic with time) to complex (for which the decay 
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Fig. 2. Plots of (a) -lm(a)/a (b) \lm(a)/Ke(a)\/q of the least stable 
axisymmetric mode against ak for various temperatures. The arrows show 
the direction of increasing T, where T/T c = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. 



of the mode is oscillatory) . 
at k = k± where 



For the least stable mode, the crossover occurs 



k: 



-(■ 



i + Vi + r 



(6) 



where q = a /(I 



a 



The parameter q increases rapidly with temperature 



and therefore k* also increases rapidly with temperature. This effect is illus- 
trated in Fig|2^a) in which, for clarity of scale, we plot —Xm{a)/a against 
(ak) of the least stable axisymmetric mode at various temperatures. The 
arrow shows the direction of increasing temperature. For k > k*, the ex- 
pression in the square root of (JSJ is positive so the normalised growth rate 
is independent of temperature and is given by Im(a)/a = $7 + v s k 2 . The 
mutual friction parameter a increases rapidly with temperature so the ax- 
isymmetric modes will decay much faster at higher temperatures than at 
lower temperatures. In Fig^b) we plot the ratio of the decay rate of the 
wave to its angular frequency, \Im(a) /1Ze(a)\ against (ak) We see that for 
all the modes \lm(a)/TZe(a)\ — > q w 0.51 at large (ak). In Fig|2fb) we 
plot \Xm(a)/1Ze(cr)\/q against (ak) of the least stable axisymmetric mode at 
various temperatures. Again we find \Im(a) /1Ze(a)\ — > q at large (ak) for 
all temperatures and that this limiting value is obtained at larger values of 
(ak) as the temperature is increased. 

In order to consider the non-axisymmetric modes we must solve the cou- 
pled equations (!2IMf) and further details of the solution method can be found 
in Henderson & Barenghi ^l In Fig. |3lwe plot —2m(a) and \2m(a) /TZe(a)\ 
against ak of the least stable computed modes for m = 0,2, 4, 7, 10 using the 
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Fig. 3. Plots of (a) the growth rate, —Im(a) (b) \2m(a) /TZe(a)\ against ak 
for the m = 0, 2, 4, 7, 10 mode. The arrow shows the direction of increasing 

in. 



same parameters as for Fig. ^ The arrow shows the direction of increasing 
m. We do not plot the m = 1 for clarity; it suffices to say that the least 
stable m = 1 mode is bound by the least stable m = mode. 

It can be seen that the modes become less stable as m is increased which 
is similar to what we found at lower temperaturea^. We can also see that 
at large (ak), all the modes are such that \2m(a)/TZe(a)\ — ► q ~ 0.51 and 
that this limit is approached from below for larger values of m. What this 
means is that oscillatory decay is more pronounced for modes with larger 
azimuthal wavenumbers, particularly at small to moderate values of the axial 
wavenumber. As for the azimuthal mode, we find that the decay rates of the 
azimuthal modes, 2m(a) are proportional to a, so at higher temperatures, 
for which a is small, the mode will decay more rapidly in time. 



3. APPLICATION 

Our work has applications to the recent discovery by Finne et aP* of 
an intrinsic condition for the onset of quantum turbulence in 3 He-B, namely 
q > 1.3. Finne et al interpreted their result in terms of Kelvin waves - helical 
perturbations of the position of a vortex core away from the unperturbed 
straight shape. In the case of an isolated vortex line, it was predicted »» 
that q = 1 is the crossover between Kelvin waves which propagate (q < 1) 
and Kelvin waves which are damped (q > 1). If q < 1 the propagating 
Kelvin waves grow in amplitude exponentially with time, driven by the local 
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difference between the normal fluid velocity and the velocity of the vortex 
line. 

If this condition is valid for a large number of vortices (which have more 
degrees of freedom that a single vortex line), then one expects that, when the 
waves' amplitude becomes of the order of the average intervortex spacing, 
the vortices reconnect which each other, and form a turbulent tangle. 

Our analysis shows that, for both axisymmetric and non-axisymmetric 
perturbations of a vortex lattice, the ratio of imaginary and real part of the 
complex growth rate is approximately equal to the parameter q identified 
by Finne et al in the large axial wavenumber limit, for all temperatures. 
The minimum permitted value of k will be governed by the height, h of the 
apparatus by the relation fc m i n = 2ir/h, so provided that the aspect ratio is 
small enough (h/a <C 1) we find agreement with the argument of Finne et al. 
Therefore if q > 1, perturbations shrink in amplitude before they can rotate 
a full cycle (overdamping) ; this reduces the volume of space which is swept 
by the vortex lines in their motion, which reduces the probability of making 
reconnections with neighbouring vortices, and without reconnections there 
is no turbulence. Vice versa, if q < 1 the perturbations propagate, which 
favours vortex reconnections, hence turbulence. 

Finally, the exact dispersion relation which we have found (12I3|) can be 
used to study with precision the torsional oscillations of on a rotating vortex 
lattice for any height and radius of experimental apparatus, which is a topic 
of current experimental intereslP. 
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